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(t) 





\ 


ct 


m 


Generalized displacement coordinate associated with the nth mode 
of vibration of body X 

Mode-dependent parameter used to define body X center-of-mass 
(see equation 48) 

Force associated with the mass of body X and its combined 
centripetal and Coriolis acceleration 

Control torque applied to symmetric wheel m 


E n x Mode-dependent parameter used to define body X inertia tensor 

(see equation 52) 

F m n x Mode-dependent parameter (see equation 70) 

Resultant force of constraint acting on body k through hinge point 
k-1 

Resultant force acting on body k due to all mechanisms existing 
between bodies J(k) and k at hinge point k-1 

Resultant of all external forces applied to body k 



Pseudo-inertia tensor of body X with respect to the origin of nest 
k-1 and the hinge point i-1 

Linear momentum of body X relative to the inertial origin 

Linear momentum of element i of body X relative to the inertial 
origin 

Inertial angular momentum of symmetric wheel m 


H m Relative angular momentum of symmetric wheel m 

H m Angular momentum of symmetric wheel m relative to the 

undeformed state of body MO(m) 

H n x Mode-dependent parameter (see equation 68) 

I w Inertia tensor of symmetric wheel m about its center-of-mass 
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V 



J(X) 



M 

MO(m) 




N 

N m 

N r 



Body label of the body to which body X is attached at hinge 
point X-l 

Mode-dependent jparameter (see equation 72) 

Inertial angular momentum of body X about the inertial origin 

Inertial angular momentum of the element i of body X about 
its own center-of-mass 

Inertial angular momentum of body X about its center-of-mass 

Total number of symmetric wheels 

Body label of the gyrostat in which symmetric wheel m is 
embedded 

Total mass of body X 
Mass of element i of body X 

Total number of rigid bodies, flexible bodies, and point masses 
Total number of flexible body modes (all flexible bodies) 

Total number of rigid and flexible bodies 
Free coordinate vector m 

Mode-dependent parameter used to define relative angular 
momentum of body X due to deformation (see equation 56) 


°r. x Position vector from hinge point X-l to the undeformed center-of- 

mass position of the element i of body X 

S D Set of all rigid-body, body labels 

S„ Set of all flexible-body, body labels 

r 

S T Set of all point-mass body labels 

S k j Set of all body labels of those bodies outboard of hinge point k-1 

relative to body 1 ; body labels of those bodies in nest k-1 

S k x t Set of all body labels associated with those bodies lying on the 

topological path from hinge point k-1 to the center-of-mass of 
body X 


vi 



SP 


ou 



t 


Tensor operator which maps vectors into skew symmetric tensors 
of rank 2 , dyads 

Position vector from hinge point A-l of body A to the center-of- 
mass of body A 

Position vector from hinge point A-l to the undeformed center-of- 
mass position of body A 

Position vector from inertial origin to hinge point 0 of body 1 ; that is, 
to the center-of-mass of body 1 
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yi,\ 

r k ,x 

5 m,„ 
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i,\ 


-> 

Aa . 


A<1> 


O 







AH 


m 


Position vector from hinge point J(X)-1 to hinge point X-l 

Position vector from hinge point J(X)-1 to the undeformed position 
of hinge point X-l 

Position vector from inertial origin to center-of-mass of body X 

Position vector from hinge point k to the center-of-mass of body X 
Skew symmetric tensor form of the vector 7 * 

Kronecker delta function 

Position vector from the deformed center-of-mass position of 
body X to the deformed center-of-mass position of the element i 

Position vector from the undeformed to the deformed center-of- 
mass position of body X 

Position vector from the undeformed to the deformed position of 
hinge point X-l 

Change in the inertia tensor of body X from the undeformed to the 
deformed state 

Angular rate of a reference frame fixed in body J(X) at hinge point 
X-l relative to a reference frame fixed at hinge point J(X)-1 in 
body J(X) 

Angular momentum of the symmetric wheel m associated with the 
deformation-induced rotation of its spin axis relative to the unde- 
formed state of body MO(m) 


vii 
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Angular velocity of a reference frame fixed in the despun wheel m 
relative to the undeformed state of body MO(m) 


Angular momentum of the flexible body X relative to its center-of- 
mass, due to the elastic deformation of all elements i of body X 

Position vector from the^yndeformed to the deformed center-of- 
mass position of the element i of body X 

Damping ratio 

Relative rotation vector of the element i of body X from the unde- 
formed to the deformed state 

Displacement about or along free coordinate vector 

Inertia tensor of body X about its center-of-mass 

Inertia tensor of element i of body X about its own center-of-mass 

Inertia tensor of element i of body X about its own center-of-mass 
in the undeformed state of body X 

Inertia tensor of element i of body X about the deformed 
center-of-mass position of body X 

Inertia tensor of body X in the undeformed state relative to the 
undeformed center-of-mass position 

Resultant torque acting on body k due to external causes 

Resultant torques of constraint acting on body k through hinge 
point k- 1 

Resultant torque of constraint acting on symmetric wheel m 


Resultant torque acting on body k due to all mechanisms existing 
between bodies J(k) and k at hinge point k-1 




The nth normal mode displacement vector for the element i of 
body X 

The nth normal mode rotation vector for the element i of body X 



<3 k Angular velocity of the body X fixed coordinate frame relative to 

the inertially fixed coordinate frame 

°w. Angular rate of a reference frame fixed in body X at hinge point 

X-l relative to a reference frame fixed in body J(X) at hinge point 
X-l 

Angular velocity of the body X fixed coordinate frame relative to 
the body J(X) fixed coordinate frame 

co w Inertial angular velocity of symmetric wheel m 


co 
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Angular velocity of the symmetric wheel m relative to the unde- 
formed state of body MO(m) 
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Relative angular velocity of symmetric wheel m 

Natural frequency of the nth normal mode of vibration for body X 

Column matrix 
Square matrix 
Rectangular matrix 

Vector it 

3X 1 matrix of components of vector it relative to body X fixed 
coordinates 


it Time derivative of vector ft relative to inertially fixed reference 

frame 

O 

it Time derivative of vector it relative to the local reference frame when 

no confusion exists as to which local reference frame the differen- 
tiation is with respect to 

1 Unit dyad 
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k-1 


Summation over all indices i contained in the set S k t 
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i@ k-1 


Summation over all indices i of vectors defined at hinge point k-1 



Sum over all elements of the body X 


£ 

m: 

MO(m)=\ 

£ 


Sum over all wheels imbedded in body X 


Sum over all modes to be used for the description of flexible body 
X elastic deformation 

Equivalence sign 
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A VECTOR-DYADIC DEVELOPMENT OF THE 
EQUATIONS OF MOTION FOR N-COUPLED 
FLEXIBLE BODIES AND POINT MASSES 


Harold P. Frisch 

Goddard Space Flight Center 


INTRODUCTION 

The development of a general purpose simulation model to be used to study the attitude 
dynamics of a general class of multibody spacecraft having one or more flexible members 
presents a number of formidable problems. The model must be computationally efficient, 
yield an accurate representation of the spacecraft's attitude dynamics, and provide suffi- 
cient generality to be applicable to a broad class of spacecraft. 

The burden of writing a computationally efficient program rests both with the analyst and 
the programmer. The analyst must write the equations in such a form that redundant com- 
putations and multiplications by zero can be readily avoided. The programmer must use 
computationally efficient programming techniques. If the programmer does not include 
the logic to avoid redundant and meaningless computation, and makes excessive unneces- 
sary use of multidimensioned arrays, the speed of computation can easily be increased by 
more than an order of magnitude above what it should be. 

The burden of defining a representative model of a complex spacecraft rests solely with the 
analyst. An overdefined model will yield accurate answers but will be useless from the 
standpoint of computational speed; an underdefined model will be highly efficient compu- 
tationally but will yield erroneous results. 

In reference 1, Velman’s nested body approach has been used to derive the equations of 
motion of a topological tree of rigid bodies, point masses, and symmetric momentum 
wheels. These equations have been programmed and form the basis for the general purpose 
digital computer program, N-BOD. Too often, during the course of a particular study, the 
question has arisen as to the effects of body flexibility and its resultant perturbing effects 
on the response characteristics of the composite system. The quest for an answer to that 
question has motivated this study. It builds directly upon the nested body approach 
presented in reference 1 . Excellent reviews of other methods of attacking the same general 
problem are provided by Likens in references 2 and 3. 

The aim of this flexible body analysis is to derive a final set of coupled body equations 
which separate effects that are analogous to those found by the coupled rigid body analysis 



from those effects which are unique to the coupled flexible body analysis. Achievement 
of this goal allows flexibility to be treated as simply an add-on effect. 

The equations are first given as a set of simultaneous vector-dyadic differential equations 
and are then reduced, by a method which eliminates unknown forces and torques of con- 
straint, to a set of simultaneous scalar differential equations. The equations are cast into 
a partitioned matrix format which has certain computationally desirable symmetry proper- 
ties. The effects of body flexibility are introduced by adding a few simple terms onto the 
coupled rigid body equations of motion, along with a generalized displacement coordinate 
equation for each elastic mode of vibration. To describe a body’s flexibility, clamped-free 
normal modes and frequencies of vibration are assumed obtainable for each flexible body. 
These may be obtained by closed-form continuum analysis, finite element analysis, or an 
educated guess. The mode shapes are not used per se within the final set of coupled body 
equations. They are used by a preprocessor to obtain resultant mode-dependent parame- 
ters for the stream of input data. 

Frequently the analyst must decide whether it is better to model a flexible appendage as a 
spring-connected rigid body or as a flexible body having several significant natural modes 
of vibration. To make an intelligent decision, there must be some means of comparing the 
two approaches. 

Because the equations derived are valid for either a rigid body or a flexible body model of 
any appendage, it is a routine process to step through the equations first with a rigid body 
model and then with a flexible body model of an arbitrary appendage. As both models 
simulate the same gyroscopic effects, equivalence relations can be established. These 
relations can then be reduced to a set of three fundamental equivalence relations, which 
provide the guide for determining the best modeling approach for any particular problem. 


BASIC SYSTEM 


The basic system to be studied is identical to that defined in reference 1 , with the excep- 
tion that now one or more of the bodies may be flexible. The theoretical development 
and notation used here is a natural extension of the work presented in reference 1 , and, 
therefore, the large volume of supportive material presented in that document will not be 
repeated. 

Consider the pair of contiguous bodies shown in figure 1. 


As in reference 1, the following definitions are made: 

Body 1 = principal body of the system of coupled bodies, 


Body J(X) = body contiguous to and, relative to body 1 , inboard of body X, 

Hinge point X-l = point of connection between bodies J(X) and X, 


a. 


= position vector from hinge point X-l to the center-of-mass of body X, 
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BODY X 



Figure 1. Notation for a pair of contiguous bodies. 

= position vector from hinge point J(A)-1 to hinge point A-l, 

<f> = inertia tensor of body A plus that of all imbedded despun 

symmetric wheels about the composite body center-of-mass, 

m, = total mass of the body A, 

A* 

7J = position vector from the inertial origin to the center-of-mass of 

body A, 

= angular velocity of a coordinate frame fixed at hinge point A-l 

A. 

in body A relative to a coordinate frame fixed at hinge point 
J(A)-1 in body J(A), and 

co^ = angular velocity of a coordinate frame fixed at hinge point A-l 

in body A relative to a coordinate frame fixed inertially. 
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The present study permits any one or all of the coupled bodies to be flexible. To maintain 
the required generality, it is assumed that a finite element model of each flexible body of 
the system can be defined. The elements into which each flexible body is decomposed 
are assumed to have both mass and rotational inertia characteristics. Of course, as the 
number of elements used increases, the rotational inertia characteristics of the elements 
approach zero. 

To establish the required notation, consider figure 2 for the arbitrary flexible body X. 



Figure 2. Notation for flexible body X. 


Let: 



= position vector from 
position of body X, 

Aa x 

= position vector from 
position of body X, 

°l x 

1 ,\ 

= position vector from 
position of element i 


hinge point X-l to the undeformed center-of-mass 


the undeformed to the deformed center-of-mass 


hinge point X-l to the undeformed center-of-mass 
of body X, 
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m 




<J> 

i,i,\ 

° 4 > 

UA 


$ 






= mass of element i of body X, 

= inertia tensor of element i of body X about its own center-of-mass, 

= inertia tensor of element i of body X about its own center-of-mass in 
the undeformed state of body X, 

= inertia tensor of the element i of body X about the deformed center-of- 
mass position of body X, 

= position vector from the undeformed to the deformed center-of-mass 
position of the element i of body X, 


5. 


i,\ 


t 


A 


= position vector from the deformed center-of-mass position of body X 
to the deformed center-of-mass position of element i, 

= relative rotation vector of the element i of body X from the undeformed 
to the deformed state, 


MO(m) = body label of the body in which symmetric wheel m is imbedded, 


H 


w 


= relative angular momentum of the symmetric wheel m imbedded in 
body MO(m), 

= inertia tensor of the symmetric wheel m about its own center-of-mass, 


% 

A 

X 

A3> 


position vector from hinge point J(X)-1 to the undeformed position of 
hinge point X- 1 , 

position vector from the undeformed to the deformed position of hinge 
point X- 1 , 

inertia tensor of body X in the undeformed state relative to the unde- 
formed center-of-mass position, 

change in the inertia tensor of body X from the undeformed to the 
deformed state, 


co. 


O 

A0, 


angular rate of a reference frame fixed in body X at hinge point X-l 
relative to a reference frame fixed in body J(X) at hinge point X-l 
(relative rate minus the effects of body J(X) elastic deformation), 

angular rate of a reference frame fixed in body J(X) at hinge point X-l 
relative to a reference frame fixed at hinge point J(X)-1 in body J(X) 
(the effect of elastic deformation of body J(X)), 
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o~ 

H m = angular momentum of the symmetric wheel m relative to the undeformed 
state of body MO(m), 

AH m = angular momentum of the symmetric wheel m associated with the 

deformation-induced rotation of its spin axis relative to the undeformed 
state of body MO(m), 

°5 W = angular velocity of the symmetric wheel m relative to the undeformed 

m state of body MO(m), and 

O 

A0 W = angular velocity of a reference frame fixed in the despun wheel m 
,T! relative to the undeformed state of body MO(m) (the effect of elastic 

deformation of body MO(m)). 

Making use of these definitions, summing over all elements i, and deleting terms which are 
quadratic in deformation, the following equations are obtained: 



= °a x + Aot x 

(1) 



(2) 



(3) 

— > 

O 

= w, + A6 

A A. 

(4) 

— > 

—y — > 



= °H + AH 

m m 

(5) 

“w 

m 

o 

(6) 

m x 

= X 

(7) 


ie\ 


o'* 

a x 

i V 

= — / m. . r. . 

m /_j 'A i,x 

X ieX 

(8) 


II 

x 3 |~ 

J3 

V 

V 

(9) 
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( 10 ) 


5. . = °r. . - °a . + e. x - Aa 

i 9 \ i,\ \ a. 


$ 

i,i,\ 


<j> 


~ ) ' m i.X S i,\ 
ie\ 

+ Ku x? m) 


= %,x + °( e! > 


= $ M,x + m i,x V^iA 1 - 6 ^ 8 ^ 


= XiA +A ^,iA + °( e2 ) 


(ID 


( 12 ) 


(13) 


(14) 




A3>. 


= ^2 ^X.i.X 
ie\ 

[(”lx -•«,)• 05) 

„,.X - ”.,x [ ( \x - ' ^X ) ’ ' ?« 1 - (V - ”“x) ? ,.x] 

+ '"'- k [ f,,x 4 (t.x - °“x ) 1 -?,.x (f tt -"“» )] 


E 


<!» 

X,i,X 


A<I> 


ie\ 


ieX 


A<3> 


\,i,X 


AH m = V * 


A0 


w 


(16) 

(17) 

(18) 
(19) 


m m 


where 


and 


^ ^ = sum over all elements i of body X, 

ie\ 

1 = unit dyad, 

0(e 2 ) = terms quadratic in deformation. 
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MOMENTUM EQUATIONS FOR A FLEXIBLE BODY 


The basic definitions for linear and angular momentum found in any textbook on rigid 
body dynamics may be applied to obtain the following relations: 

• Cjj x = linear momentum of flexible body X relative to the inertial origin 



£ 

ie\ 




X ^I,\ 


( 20 ) 


• L, = angular momentum of flexible body X relative to the inertial origin 




J m * 


H 


MO(m)=\ 


where 


( 21 ) 


= sum over all wheels imbedded in body X, 

m: 

MO(m)=\ 

( 22 ) 
( 23 ) 

and the closed dot implies an inertial derivative while the open dot implies dif- 
ferentiation relative to the local body X fixed reference frame. 

Direct substitution and application of previously made definitions yields: 




= w- + e. 

\ i,\ 


,X 


— X 6. - +5. - 
x i,\ 




= 4\ 


“x + n V 


I,x 




+ AL K\ + 




MO(m)=X 


( 24 ) 


where 

At, x x = angular momentum of the flexible body X relative to its deformed 

center-of-mass position due to the elastic deformation of all elements i 
of body X 
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(25) 


AL x,x = 


Y Ki,x‘\A +m iA^ x v| • 

ie\ L J 


H = angular momentum of symmetric wheel m, relative to the inertial origin 


H m =I W *^ + H m 

m 


(26) 


Straightforward differentiation yields the rates of change of momentum relative to the 
inertial reference: 


G i,x =m \\\ 


(27) 


L i,x = *x' V + m x7i A X7 I>x 

+ A V + S [h + co. X H 
L m \ m 
m : 

MO(m)=\ 

+ *x • «x + "x x ^x • “x + A \,^ 


(28) 


and 


9.-V (L -S^hJ 

m ' m / 


(29) 


EQUATIONS OF VIBRATION FOR A FLEXIBLE BODY 

A set of equations which define the vibrational motion of an arbitrary flexible body rela- 
tive to its undeformed state is required for the complete definition of all terms contained 
in the momentum equations (27, 28, and 29). It is not essential, but it is most convenient, 
to work toward a set of equations which define deformation in terms of the flexible body’s 
normal modes of vibration. 
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At the hinge point A-l of the flexible body A, a body-fixed reference frame is defined. 
Body A is assumed to be clamped in both translation and rotation at the origin of this 
body-fixed reference frame; all elastic deformation is measured relative to it. While the 
position and orientation of this body A fixed reference is independent of deformation, 
it may have as many as three degrees of rotational freedom relative to its contiguous 
body, body J(A). 

Each element i of the finite element model of the flexible body A has, at most, six degrees 
of freedom: three translational and three rotational. In addition, the elements are each 
subject to viscoelastic loads and conditions of constraint. 

Let 



linear momentum of element i of body A relative to an inertial reference 



( 30 ) 


and 


L. . . = angular momentum of the element i of body A about its own center-of-mass 

1,1, A 

relative to an inertially fixed reference frame (small angle deformation 
assumed) 



( 31 ) 


The equations of motion for the elements of body A are obtained by defining the holonomic 
conditions of constraint which limit the number of allowable degrees of freedom for each 
element and by equating the time rate of change of momentum for each element to its 
viscoelastic loading. 

In order to separate relative-translational and relative-rotational deformation, partitioned 
matrix rotation is used to write 
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' 2 ,\ 


> = < 


1 ,\ 


2 ,\ 


( 32 ) 


The holonomic conditions of constraint can, in general, be written in the form 


LgJ 




( 33 ) 


where the number of rows of the rectangular matrix LgJ is equal to the number of 
constraint conditions and hence the number of degrees of freedom to be suppressed. 


The inertial time rate of change of momentum can, in general, be written in the form 


— > 
G 

d 1 

ml 


t + ? 1 


. = — 




-> 

dt 



o 

L 




co + 6 


( 34 ) 


where row i of the upper partition is row i of the lower partition is L { . x , and 

d l /dt is the time derivative relative to the inertial reference. 


The viscoelastic conditions existing between contiguous elements can be expressed in terms 
of resultant viscoelastic loads. In general, this will take the form 



( 35 ) 


where row i of the upper partition defines the resultant viscoelastic force acting on the 
element i, and row i of the lower partition defines the resultant viscoelastic torque acting 
on the element i. 


I 


ii 


By a direct application of the' preceding, the general form for the equation of relative 
elastic deformation is found to be as follows: 


d 1 

ml 

0 


’ -V ^ 
7+5 


c n 

Cn' 


° 

€ 


k n 

k 12 


-> 

e 


0 

— 





+ 





+ 





> - . 


dt 

0 

$ 


w + e 


_ C 21 

C 22_ 


— >■ 
e 


_ k 2i 

k 22_ 


Q 


0 


subject to the constraint equation 



( 36 ) 


(33) 


In order to put the deformation equation into a form which is amenable to an eigenvalue 
solution, it may be rewritten as 


ml 

0 


OO • 
-> 

€ 


c n 

1 

<N 

u 


o • 
-V 

e 


k n 

1 

<N 

& 


e 


f t 





+ 




o 

+ 





. = * 





— y 





—>■ 





-> 


f r 

_0 



d 


n 

to 

^22 _ 


e 


_ k 2 1 

k 22- 
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where, on the right-hand side of equation 37, 


pT I 
p'R I 


ml 0 
0 <t> 


J O O 


d 1 

ml 0 


-* ^ | 
7 + 6 

dt 



O 1 


o 


co + 6 j 


(37) 


(38) 


and the constraint equation (33) 



remains unchanged. 

An eigenvalue solution is obtained by assuming a solution of the form 


(33) 



to the load-free equation 


ml 0 


O O ■ 
€ 


— i 

JT 

J* 4 

to 

1 


~> 

e 


0 




+ 




. = 




% 




-> 



_0 <t> 


e 


Al k 22_ 


0 


0 


(39) 


( 40 ) 
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subject to the conditions of constraint 



Let 

o> n = natural frequency of the nth normal mode of vibration of body X, 

A 

> = nth orthogonal mode of vibration of body X, 
and impose the following normalization conditions: 



( 33 ) 


Eh 

ieA 


~>T 

A m,i,A 


->T , ->R 

ip 1 . % . 

^ n,i, A ^ m,i,\ 


C*u,x • ?„ R ,u)] 


= rru 5 

A m ,n 


where 

x = rh/z normal mode displacement vector for the element i of body A, 
w R . , = rh/z normal mode rotation vector for the element i of body X, 

r n,i,A 

and 

5 = Kronecker delta function. 

m , n 


(41) 


It should be noted that, in the normalization equation, as the number of elements used 
increases and their respective mass properties decrease in magnitude, the contributions due 
to their respective rotational inertia effects approach zero. Hence, for all m and n, 


Lira 
N. . 

i,\ 


E 


Y m ,1, A 



A 


-f R 
^ n,i,A 


) =0 


( 42 ) 


where 

N. x = total number of elements i of body X. 

To obtain a modal solution to the deformation equation, a solution of the form 

n,X 
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(44) 





n,A 


is assumed for the motion of the element i of body X, where 

a n x (t) = generalized displacement coordinate associated with the nth mode of 
vibration of body X 

and 


£ 


= sum over all modes to be used for the description of flexible body X 
elastic deformation. 


Substitution of the assumed solution into the deformation equation (37), along with the 
standard application of orthogonality conditions, readily yields the uncoupled equations 
of motion for the generalized displacement coordinate a n ^(t). That is, 


m. 


[a„,x( t ) + 2?co nA a nX (t) + co2 


X a n,A v 


(t) = 


£ 


ieX 


->T 
^ n,i 




n,i,X 1 i,X 


(45) 


where proportional damping has been assumed and f is the damping ratio. 


Two points should be especially noted. First, each eigenvector satisfies the constraint 
equation (33); that is, 



<A 

^n,X 


| 0 | for n = 1 , 2, . . . 


(46) 


Secondly, it is not necessary to actually perform the finite element analysis outlined for 
every problem. The specific method used to obtain mode shapes is immaterial. The 
boundary conditions and normalization procedures, however, must be adhered to as 
follows: 


• Boundary conditions 

The flexible body must be clamped in both translation and rotation at one and 
at only one point; this point is the origin of the body-fixed reference frame at 
the hinge point. Reference frames fixed at other points are assumed free to 
both translate and rotate relative to it. 
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Normalization of modes 

A somewhat unconventional normalization procedure was found useful in 
uncoupling the generalized displacement coordinate equations. Six parameter 



mode shapes are assumed available: three translation components of for 

the element i and three rotation components of i p for the minicoordinate 

frame fixed in the element i. They are normalized as follows: 


E[ 

ieX 




I?* 

^ n,i,\ ^m,i,X 


(X.,x •?,«)] * 


m.5 
X m,n 


(41) 


The problem is essentially solved at this point. A tedious exercise in vector- tensor opera- 
tions must now be performed to put the equations together into a computationally effi- 
cient form. The reader is referred to reference 1 for supportive material and to references 
4 and 5 for proofs of the vector-tensor identities used and listed in the Appendix. 

There is no point in carrying out summations over all elements i at every integration step 
All summations of this type may be replaced by resultant mode-dependent parameters as 
follows: 

• Center-of-mass vector 


“x = a x + 


where 


~r E -UA.-V £ w> 

X ieX n,X 


~T ^ -*T 

An ’ x = m. 2-/ ”*•**“•*•» 


ieX 


Inertia tensor 


®x* £ [x.u + i *x.u] -X + £ E .,xVx« 

ieX n,X 


where 

D 


W"H.x [Pu-^x) •?L.x1-(\x-^)?L,x]. 


(47) 


(48) 


(49) 


(50) 


D n,x ~ yi D n,i,x ’ 

ieX 

and using the post superscript T to imply “transpose” 


E 


n,X 


= D 


n,X 


+ D 


T 

n,X * 


(51) 


(52) 
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• Angular momentum due to deformation 


ie\ n,\ 


where 


i XT' <r> 
B ”A-^ 2 ^ -‘A 

* ie\ 

Si.X i,i,\ ^n,i,X ■ 


x X *»,i,* 


ieX 


^n,X = ^n,X + m X (®n,X-°“x X ^n,x) > 

and terms quadratic in deformation have been deleted. 
• Generalized force 


( 53 ) 


(54) 


(55) 

(56) 



ieX 



+ </> 


n,i,\ 



(57) 


From the partition matrix equation (38), 

00 / •* 

6 i,X “ \^I,A + \\ 

and 




O T 

. 

1,1, A 


o o 



d 1 

dt 



(58) 


(59) 


It is convenient, in the process of removing all summations over the elements i, to separate 
and collect terms which can be associated with the various components of acceleration. 
Accordingly, write 


yi,x + 8 i,x-\x- 


«x + 


o-» < 

“x + 


r i,x- “x + e i,x 


(60) 


and carry out the prescribed differentiations in equations 58 and 59 to obtain 


% ■ -U [?.,X 4 *' *• 4 * - V\x) 


( 61 ) 
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and 


R 

i,X 



( 62 ) 


Consistent application of the initial assumption that terms quadratic in deformation are 
negligible yields the following components of the generalized force. They are obtained 
from a direct substitution of equations 61 and 62 into equation 57. 


• Generalized force component which is due to the linear acceleration of the 
undeformed center-of-mass position of the body X 


’ 2 m i,x + = (\x-“x + °“x) ‘V (63) 

ie\ 

Generalized force component which is due to the inertial angular acceleration 
of the reference frame fixed at the origin (hinge point X-l of body X) 


-£ 

ieX 

— >■ 

= - w x 

where (3 n K is defined by equation 56. 

• Generalized force component which is due to the inertial angular velocity of the 
reference frame fixed at the origin (hinge point X-l of body X) 



- ^2 j m i,\ [“x x K x (\x ~ °“x + %x)) + 2 *x x V] * *n,u 

ie\ | 

+ fix-.., • (5 X + tA~] x 

X 1,1 , X \ \ i,X ) \ ^n,i,\ 


= OK 


( H n,X + ^2 F "*.n.X a m.X (t) ) 


2w x 


^m.X® K m,n,X 


m,\ 


(65) 


where, by making use of equations 43 and 44 


?,,x= 5Z ' "' :?T 

n,X 




(43) 
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(44) 



VxW* 


,-tR 

n,i,X 


and a manipulation of vector-tensor identities, the following relationships are 
established: 


" S I ["* X ( X ( ***•* " °“ x ))] ’ + ^ X Xw * 

ieX I 


= W X* H n,x' W X 


( 66 ) 


where 


Vx- E 


and 


ie\ 


H x = D x + J x 

n,A. n,\ n,A 


(67) 


( 68 ) 


b. - Y' ■»„ [S K x (S,x?J] -C,x - 3 • V 


La 

ie\ 

where 


m,\ 


F m,n,X a m,xW ’ W X 


_ X ^ r ->T ->T 1 ->"T ~>T "I 

m,n,\ ~~ m i,X ~ ^ m,i,\ ^n,i,\ J 

ieX 

- ^2 j 2 m i,x ( -X X «i,X ) * *n,i,X + (X X Xu * r u ] * ?Ii.> 

ie\ l 

= - 2 c °X* a m,X^ K m,n,X 


(69) 


(70) 


(71) 


m,\ 


where 
— >• 

^m,n,X 


- E j m u ?«,u x ? »,u + 7 (°*u,x-?«.x) x ?»V 

ieX | 


' (72) 
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A direct summation of the above components yields the desired form of the generalized 
force for the nth mode. That is, 

Y |/n,i,X * % + *n,i,X * %] = - m X (^I,X “ K + °«X ) * \,X ~ G n,X * “x 


ie\ 


+ CO^ 


H n.X + F m,«,X a mA (t) 

m,\ 


■3,- A - 2_, (73) 

m ,X 




These definitions, established by equations 47 through 73, are used to rewrite the rate of 
change-of-momentum equations (27, 28, and 29) as 


G i,x = m x7i,x> 

*1,X = <V W X + m x 7i,x X 7 ia + G n,X a n,X^) 


(74) 


n,\ 


m : ' f 


+ aJ, X H_ 1 + <J\ • co^ + X | <3\ • co^ + AL^ 


m * m ' ~ \ X w \ “X,X 


) , (75) 


m : 
MO(m)=X 


K = V -4 + H m +^X (i w -^ + H m ) , 

m \ m / 

and the generalized displacement coordinate equation (44) as 
m x [ a n,xW + 2 f W n,X a n,xW + W n.X a „,^ ] 

= - m X \7i,x-“x + "x j •An.x-Qn.x’"; 

/ 

E 


(76) 


+ cu • H„ , + 


\ J ' *n,X ‘ / ^ F m,n,\ a m,\^ I ' 

m A / 


-2w. 


X ^ v ^m,X^ K m,n,X 
m,\ 


(77) 


In summary, the resultant mode-dependent parameters used in the analysis are: 


n,X 


nu 


Y m i,X^n,i,X 


(48) 


ieX 
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(50) 


D n,i,X = m i ; 


n,X 


"n,\ 


B 


n,A. 


n,\ 


n,\ 


n,\ 


n,\ 


m,n,\ 


m,n,\ 


,*[( (v*.) ?ij 

E 


D . . 
n,i y A. 


= D n,X + D i,X 


“ ^ m i,X V X ?Ji,X 

X ieX 

E Xi,x -? M,x 

(S„,X-^X J n.x) 

E 


ie\ 


C n,X + m X 


X 3> 

,\ A i,i,\ 


ie\ 


D n,X + J n,X 


(51) 

(52) 

(54) 

(55) 

(56) 

(67) 

( 68 ) 

(70) 


ieX 

- E ('v^ x ^j (72) 


ie\ 


NESTED BODY EQUATIONS OF MOTION 

In reference 1 , the nested body approach was used to write the coupled rigid body equa- 
tions of motion. This same approach can be used if the bodies are flexible. The extension 
of the equations to account for body flexibility simply involves carrying along a few more 
terms. The basic form of the equations remains unchanged. 

As in reference 1 , several sets of body labels must be defined. The bodies are each given 
distinct integer labels such that, along the topological path from body 1 to any limb end, 
the body labels are of increasing numerical magnitude. Let 

S = the set of all body labels associated with rigid and flexible bodies, 

S. = the set of all body labels associated with point masses, 
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S F = the set of all body labels associated with flexible bodies, 

S k l = the set of all body labels outboard, relative to body 1 , of hinge point 

k-1 of body k (these body labels define the bodies contained in the nest 
k-1), and 

S k .i x .j = the set of all body labels associated with those bodies on the topological 
path outboard from hinge point k-1 to the center-of-mass of body X. 

Furthermore, as in reference 1 , let 

jifj = position vector from the inertial origin to the hinge point of body 1 , which, 

by definition, is its center-of-mass, and 

7 " k x = position vector from hinge point k-1 to the center-of-mass of body X, 
where 


?i,x =^Wo,x ’ (78) 

y k .i,\ = y ; i+a- ( 79 > 

ieS k-l,\-l 

i*k 

? k i = resultant force of constraint acting on body k through hinge point k-1 , 

= resultant torque of constraint acting on body k through hinge point k-1, 

0^ = resultant torque of constraint acting on the symmetric wheel m which 

prevents motion about any axis normal to its spin axis, 

? k j = resultant force acting on body k due to all mechanisms existing between 
bodies J(k) and k at hinge point k-1, 

^k i = resu lt an t torque acting on body k due to all mechanisms existing between 
bodies J(k) and k at hinge point k-1, 

CL m = resultant torque acting on the symmetric wheel m due to all mechanisms 
existing between it and the body MO(m) in which it is embedded, 

= resultant external force acting on body k, and 
0^ e > = resultant external torque acting on body k. 
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Making use of the above definitions and Newton’s fundamental laws of motion, the equa- 
tions of motion for the nest k-1 are: 


E K - E 


(80) 


AeS 


k-1 


\eS 


k-1 


2 iu fa.-**) * (*l.**&) 

tes k-l 

^ s k-l ' ’ 


(81) 


and 


— > — >■ — > 

H=** + CL 

m W m 

m 


(82) 


By direct elimination of the forces of constraint from the angular momentum equations, 
the following relations can be shown: 


• Rotation equation for nest k-1, keS R 

E (V (*•«-*) 


AeS 


k-1 


Xes k-1 


(83) 


Translation equation for nest k-1, keS R and keS L 


XeSk.j XeS k-l 




(84) 


Symmetric wheel equation for wheel m 


H m =*w + CL m 

m 


(82) 
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• Generalized displacement coordinate equation 


m X [vxW + 2fW n,X * n ,x(t) + W n,X a n,xW] 

(\x-“x + 5 “x) ‘An.x-Qn, 

' (»n,X + Z 


= -nu 




F m n x a m 

m , n , A. m 


*m # A^ ^m,n,X 


,x«>) 


U, 


(85) 


where, from equations 74, 75, and 76, 

G I,X = m x' ), i,x ’ 

Li,x" (\k-“k) XG i,x = ^•“^"'xlu^U 


(86) 


Q . a . + 

n,\ 


2 (»n, +3 x x »,) 


n,\ 


m : 

MO(m)=\ 


+ \ ’ w x + co x x K‘ W X + AL X, 


(<t>x • + ALx ( x) , 


and 


H_ = I„, • 3, +H„ +«, X I w • co, + 


*m W w \ m W X 

m 


K- s ^) 


(87) 


(88) 


Special consideration must be directed toward the vector representation of the translational 
acceleration of the center-of-mass of the body X. Since body X may be at the end of a 
chain of flexible bodies, several terms which were zero in reference 1 must now be included. 
That is, from equations 78 and 79, 

?i,x = K [^j(i) X 0i + "j(i) X X Pi) + 2 ^J(i) X + 13 i J 


i*l 


+ w x x “x + w x x 


K x «x) 


+ 2w x Xa x + Q: x • 


(89) 
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From the standpoint of optimizing computational efficiency, it is desirable to establish 
symmetry relations. Accordingly, as in reference 1 , all inertial angular acceleration vectors 
are replaced by summations of relative angular acceleration vectors. The identity 

(90) 



is therefore introduced and equation 89 rewritten as 

i*-f,*12 ,x 


ieS 


0,\-l 


ieS 


0,\-l 


^ l [ W J(i) X (^J(i) X @i) + 2C0 J(i) X ] + 03 \ X ( W \ X “x ) 


+ 2w x X a x . 


ieS 


0,\-l 

i#l 


(91) 


Substituting the various vector identities provided and recombining terms yields: 
• Rotation equation for the nest k-1 , keS R 


XeS 


k-1 


12 12 


(i) 


m n,\ 


MO(m)=X 


ieS 


0,k- 1 


i^k 


y wy: ^-i,k xm x ( e5 i x ^i,k) + ^2 ^-i^ xm Ji 

XeSk.j ^ S k-1 XeS k-l 




ieS 


k-1 


y \ ^ + y >k.i,x x m x ("i x ^ i.x) + y ^ k ->> 

\eS._ j XeS i-l XeS. 

\*s T 


X x m x /3.. 


\eS 


k-1 L 


V “x + “x X K • "x + A \,x) + y; ^XH, 

m: 

MO(m)=\ 
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?k-i,x Xm x (^« X (^0) X ^) 


+2 “j(i) x 0i 


+ to x X 


fa x «x) 


+ 2 “x x «x 




Translation equation for the nest k-1, keS R and keS L 

E m *fa + “J + S E m X^i X ?i-l,X + E 

* eS k-l ieS 0,k-l 1 XeS k-l * eS k-l 


E (' 


w j(i) x K ( i) x ^i) +2cj j ( i) x ^i 


(5„xs„) 


+ x “a + f ", +f J. 


.♦E 


Symmetric wheel equation for wheel m, X = MO (m) 


vz 


! w * «, = -«„ X I w -w„ 


*H m )* 


CL +0J 


Generalized displacement coordinate equation 


I OO 

l n,X (t)+ ^l *K»,X + E 


E ^(^iX^-l.x) * \,X + ^n,X*^il 

*s h J 


25 



= -ITU 


-m. 


[ 2 K,^Vx (t) + W a a n,X (t) ] 

yi (^ja) x (^j(o x ^i) +2 ^(i) x + ^\ x ( w \ x a \) 

»eS n . , 


n,\ 


Ol 


+ C0, 


2cX 


^n,\ + ^m,n,\ a m,A^ 


m,\ 




^ a mA (t) ^m,n,X 


(95) 


m,\ 


where the second order effect assumption has been exercised to write 

y; s i x (^.i.\ _ ^ +d ®x) = s s i x ^-i,x 


(96) 


ieS o,x-i 

in equation 95. 


o,x-i 


These equations may now be put into their final vector tensor form by the introduction of 
several more definitions. As in reference 1, let: 


SP= skew operator such that for vector v" having components (v t , v 2 , v 3 ), 
&(*) = 


0 V 3 ' V 2 


-v 3 0 V, 

L V 2 ' V 1 °. 


and 


Furthermore, define 


G k-i ( i-i =m x * ?k-i,x) ’"Ww] 

^k-l.X = '^ ? ( 7 k-l,x) • 


(97) 


(98) 


( 99 ) 


? = force associated with the centripetal plus Coriolis acceleration of the 

center-of-mass of body X relative to the inertial reference 
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(100) 


t 


x,x 


c x = m x 


23 ( I? « )< ( 3 I(1I*?) ,!3 I» X M 


ieS 


0,A-1 

i#l 


+ “x x 


("x X “x) + 2*x X “x) 


angular momentum of body X about its own center-of-mass relative to 
the inertial reference 


ft 


L x,x -<i, x * CJ x + AL x,x + 


H 


( 101 ) 


m: 

MO(m)=\ 


angular momentum of wheel m about its own center-of-mass relative to 
the inertial reference 


H - I w • co, + H 

m W Am 


(102) 


and make use of the vector-tensor identities provided in the Appendix to recognize that 

^k-i,x X - - lVi,x ‘ , (103) 

t k . 1 > x Xm x (“i X ti,x) =G k-i,i-i *"i ’ (104) 

I I * (105) 

m W W 5 

m m 

and, 

y •“> — > 

^i-l,A * ^n,A “ * ^n,\ * ( 106 ) 

Substitution and rearrangement of terms yields the following relations: 
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Rotation equation for nest k-1, keS R 


ieS 


0,k-l 

i#k 


Z*- Z 


k-1 ,i-l 


XeS^ j XtS 

^s L 


k-1 


00 . 


ieS 


k-1 


ieS t 


Z^Z 

XeS. t A.£S| j 


k-1 ,i-l 


i-1 

\^s T 


CO. 


Y Uwd + [ _rn ^ r Ki,x '^il 

^S. , L J ^S k ., L 


k-1 
\eS T 


*Z 

° o 

Y _m k r k-iA'^ 

*Z 

00 

Yj " m X r k-X,X '% 

ieS 0,k-i 

XeS k-l 

ies k-l 


i^k 
i=£ 1 

*Z 

Z V“». + Z f 

T V r k-l,x'^n,X + ^a) *n,X 

XeS k-i 

m: 

n,\ 



L MO(m)=X 

\eSp 



£[ 

Xe5 k .i L 


< V W X + W X XL X,X + 7 k- 






s[ 

Xe3 k .i *" 


i,^ e)+ ¥ 


( 107 ) 



• Translation equation for nest k-1, keS R and keS L 


X ' oo X > X ' -*• 

m x 1 ' + 2-r 1 * “x + Z Z™* 1 ’ A "> x a "> x (t) 

\e s k i XeS k-l XeS k-l 

\eS T 


\eS r 


E 


kS 0,\-l 
i#k 


Z r »-»A ‘ + Z 1 


XeS 


k-1 


XeS k-l 

i=£l 


ieS 


k-1 


Z r M A * + 2 m X 1 


\eS. , 
l-l 


\eS. 


i-1 


2 2 ? 

XeS k-l XcS k-l 


(e) 

X 


• Symmetric wheel equation for wheel m, X = MO (m) 


* W i + *W m ’ W w ■ W \ XH m + CL m + 
m mm m 


ieS 


0,\-l 


• Generalized displacement coordinate equation for mode n of body X, XeS r 


E (- 

* eS 0,\-l 


\ i 

n,xj ' °°i 


" m X r i-l,X • A n,X + Qn.X • + m X A n,X * 


+ m x2 n ,x( t ) + m x A nX 


_ o o 

E* 


ieS 


o,\-i 

i*l 


= -m, 


-nu 


a n,X (t >] 

Z ( 5j d) X (^(0 X ) + 2 “j(0 X ) + “x x (“x X °\) 


ieS o,\-i 

L i=£l 


+ CO, 


• H n,X + Z F ">."A a m,\W 

m , A. 


“x" 2 “x 


E s 


a m x(t)K, 


m,\ v ' m,n,\ 


m,\ 


(108) 


(109) 


n,\ 
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PARTITIONED MATRIX FORM OF EQUATIONS OF MOTION 

The preceding equations may now be put into a compact partitioned matrix format. To 
do this, the number of bodies, wheels, and flexible modes must be known. Let 

N = total number of rigid bodies, flexible bodies, and point masses, 

M = total number of symmetric wheels, 

N m = total number of flexible body modes, and 

N_ = total number of rigid and flexible bodies. 

K 


By inspection,* it is readily seen that the equations of motion, equations 107 through 1 10, 
may all be expressed by the single partitioned matrix equation 


X 

i c 

F 

r 


• N 

GO 




<P° 


' 0< e > ' 

T 

I C 

i s 

0 

0 


o 

w w 


7 2 




CL 





* 4 


► = 


► + < 


► + i 

► 

f t 

0 

M 1 

A 


a 


% 


0 


0 

r T 

0 

A t 

M 2 


O 0 

0 I 

k y 


. ^4 . 


F c 


p( e ) 

h. J 


( 111 ) 


where the coefficient matrix on the left-hand side of the equation is symmetric (post 
superscript T implies “transpose”) and the following definitions are applied: 

• X • co + I c • co w + Fa + f • 0 =rt 1 + 0 C + tf> (e) (H2) 


This is an (N+1)X 1 column matrix of vector equations defining the equation of 
rotation for each nest containing one or more rigid or flexible bodies, the equa- 
tion of translation for each nest containing a single point mass, and the equation 
of translation for the composite system. 


I c • co + I s 




+ CL 


(113) 


This is an MX 1 column matrix of vector equations defining the equation of 
motion for each symmetric wheel. 

F T • w + M 1 a + A • (3° = tj 3 (114) 

This is an N m X 1 column matrix of scalar equations defining the generalized 
displacement coordinate equation for each mode of vibration of each flexible 
body. 


*The reader is referred to the section “Equations of Motion (Matrix Form)” in reference 1 for supportive material- 
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(115) 


• r T • co +A T a + M 2 • < H > = j? 4 + F c + Fl e > 

This is an(N R -l)X 1 column matrix of vector equations defining the translation 
equation for each nest containing one or more rigid or flexible bodies but not 
the translation equation for the composite system. Since relative translation is, 
by definition, not permitted between contiguous rigid or flexible bodies, these 
equations do not introduce any independent state variables. Their inclusion, 
however, provides certain computationally useful symmetry relations and the 
equations which define the unknown forces of constraint acting between con- 
tiguous rigid or flexible bodies. 

The partitioned matrix equation defined above makes use of an expanded set of state varia- 
bles which are not all independent. In order to arrive at a set of simultaneous equations 
which will not lead to singular coefficient matrices, a set of linearly independent state varia- 
bles and a rectangular linear transformation matrix must be defined. 

Throughout the development, an attempt has been made to separate those effects which 
are obtained by a straightforward coupled rigid body analysis from those which must be 
added on to account for the effects of body flexibility. Accordingly, it is convenient to 
choose as the set of independent state variables those rate vectors which define the relative 
motion of contiguous bodies at hinge points and the generalized displacement coordinate 
rates which define the effects of flexible body vibration. It should be noted that a some- 
what loose definition of independence is being employed at this point as the relative rate 
vectors used coalesce as many as three degrees of relative freedom into a single independent 
state variable. 

Making use of previously made definitions, the various elements of the expanded set of 
state variables may be expressed as follows: 

• Relative angular rate between the reference frame fixed in body X 
and the reference frame fixed in flexible body J(A) 

o 

^ 

co^ - °gj>^ + A6 X XeS R (116) 

and 

o 

= XeS L (117) 

where 

o 

A0 X = Vj ^n,J(X)( t )^^,H(X-l), J(X) ’ (118) 

n,J(X) 
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^iTh(x-i) J(\) = nt ^ norrna l mode rotation vector at hinge point A-l on the 

flexible body J(A), 

and 

o o o 

"x-°«x<A + *x x 4 • (11 

• Relative translation rates 

For all point masses contiguous to either rigid or flexible bodies, the relative 
translation rate 0 

*X XeS L 

is an independent state variable. Similarly, the translation rate of hinge point 0 
relative to the inertial reference 


% 

is also an independent state variable. 

• Relative angular rate between symmetric wheel m and flexible body MO(m) 
fixed reference frames 


where 


w w w w + 

m m nu 


Aft 


w 




n,MO(m) x ^ ^n,W(m), MO(m) 9 


n, MO(m) 


( 120 ) 


( 121 ) 


<2f R ,„. , , = nth normal mode rotation vector at the attachment point 

of the wheel m in the flexible body MO(m), 


and 


co. 


w 


= co. 


w 


+ Aft 


w 


( 122 ) 


• Translational rate of hinge point X-l relative to its undeformed position in 
flexible body J(A) 


K = (123) 
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where 
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^i£h(x-i),j(x) = nt ^ n o rm al mode displacement vector at hinge point A-l in 
the flexible body J(X), and 


/3x = ^, 


(125) 


Combining the above definitions, it is possible to write the following rectangular linear 
transformation matrix equation that relates the expanded set of state variables to the 
chosen set of linearly independent state variables. That is, in partitioned matrix form, 
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Direct substitution of equation 127 into 111 yields the final form of the simultaneous 
vector-dyadic equations of motion for the total system and the equations which define 
the internal forces of constraint. That is, 
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By leaving the equations in this form it is obvious that the effects of body flexibility 
may be taken into account by simply adding a few extra terms to the coupled rigid 
body equations. In practical applications, a simplified coupled rigid body model is 
usually the starting point. At this level there are no flexible modes of vibration and 
the equations reduce to those derived in reference 1. In the next step, some of the 
bodies are taken to be flexible, but hinge points and wheel attachment points are assumed 
to be either at node points of flexible bodies or in rigid bodies. At this level the elements 
of the partitions , <p* , and are all zero. In the final step of analysis, both hinge 
points and wheel attachment points are subject to the effects of flexibility. At this level 
all partitions have nontrivial entries and must be included. Computationally, the above 
form is advantageous since simple logic flags can be created from an interrogation of the 
input data. These can then be used to determine whether various computation loops 
should or should not be executed. 


ELIMINATION OF CONSTRAINT TORQUES 

To solve the coupled vector-dyadic equations of motion, the unknown constraint torques 
must be deleted and the equations reduced to a set of simultaneous scalar equations. The 
procedure is identical to that used in reference 1 . 

At every hinge point, free coordinate vectors are defined. These define the vector direc- 
tions about or along which relative motion of contiguous bodies is allowed (free). These 
may be sequentially indexed and the relative angular velocity vectors expressed as 



(129) 
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where 


cf m = free coordinate vector defined at hinge point k-1, 

» 

0 m = scalar rate of ,body k relative to body J(k) about or along the free vector 
q m ,and 



= sum over all free coordinate indices m, defined at hinge point k-1. 


m @k-l 


Differentiation of equation 129 yields 
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Similarly, for symmetric wheels we may write 


and for point masses, 
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A partitioned rectangular matrix can be defined such that the system state vector can be 
rewritten as 
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Furthermore, at each hinge point, the free coordinate vectors are orthogonal to the con- 
straint torque, hence 
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In a manner identical to that used in reference 1 , the following set of simultaneous scalar 
equations is obtained by direct substitution of equation 134 into 128 and application of 
the orthogonality relation defined above in equation 135: 
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where the last N R -1 vector equations of the partitioned matrix equation (136) may be 
evaluated if the forces of constraint between contiguous rigid or flexible bodies must be 
known. 

Equation 136 implies that a horrendous amount of matrix multiplication must be carried 
out. From a computational standpoint, however, the matrices are sparse, and relatively 
simple logic statements can be created which will eliminate all unnecessary numerical 
multiplications by zero. 


PREDEFINED RELATIVE MOTION OF CONTIGUOUS BODIES 


In practical application, it is frequently desirable to model an appendage deployment or 
the relative motion of a scanning antenna without introducing the details of a high fre- 
quency onboard control system. The ability to do this often results in an enormous saving 
of computer simulation time, since the high frequency loops of the control system need 
not be modeled. The solution method is rather straightforward and can readily be out- 
lined in general terms. 

The simultaneous vector-dyadic equations which define the motion of the system of 
coupled bodies take on the general form 

[Z] {w} = \u\ + \ <P C } (137) 


where, from equation 1 28, 
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If one or more of the bodies of the system has its relative motion defined, the number of 
independent state variables is reduced, additional unknown constraint torques are intro- 
duced, and special solution techniques must be used. Let 

co = column matrix of all relative acceleration vectors which are given by known 

c 

functions of time and state, 

<% = column matrix of all independent relative acceleration vectors, and 

[R] = reordering matrix such that 


rss 



( 142 ) 


Direct substitution of the identity into the basic form and premultiplication by the trans- 
pose of the reordering matrix yields 


where 



( 143 ) 


i 0* l = column matrix of constraint torques acting between rigid or flexible bodies 
' ' and constraint forces acting between point masses and their contiguous 

bodies at hinge points at which relative motion is undefined, and 

) 0* l = column matrix of constraint torques acting between rigid or flexible bodies 
' ' and constraint forces acting between point masses and their contiguous 

bodies at hinge points at which relative motion is defined. 

Since all elements of the column matrix 
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are given functions of time and state, the equations of motion for the coupled body system 
are defined by 


M • (“4 * M * (<} - M • N • 


(144) 


and, if desired, the unknown constraint torques can be obtained from the equation 

kl'M [ z 4 <145) 

FLEXIBLE VERSUS RIGID BODY MODELING 

In establishing a representative model of a spacecraft with flexible appendages, the analyst 
must decide whether to model each appendage as a spring-connected rigid body or as a 
flexible body having a finite number of significant modes of vibration. 

Neither modeling technique is perfect. Both have advantages and limitations which, if 
not taken into account, can severely compromise the conclusions drawn from the simula- 
tion model. It is possible, by a detailed study of the analysis presented in this report, to 
assess the analyst’s choice of modeling technique. This is done by deriving a series of 
equivalence relations which define how the two modeling approaches describe the funda- 
mental dynamic characteristics of the elastic appendage and their dynamic coupling into 
the system equations of motion. 

Whether the analyst chooses to model a flexible appendage as a spring-connected rigid 
body or as a flexible body, certain fundamental dynamic effects exist and must be proper- 
ly accounted for in the system equations of motion. These fundamental dynamic effects 
are independent of the choice of modeling technique; hence, the terms which simulate 
them, as derived by the two modeling approaches, are equivalent. 

If the topological tree of connected bodies has a flexible body which is not at a limb end, 
terms appear in the equations of motion which cannot be simulated by a simple spring- 
connected rigid body model of that flexible body. These terms are related to the fact that 
coordinate frames at the hinge points on the flexible body are not fixed relative to each 
other. 

The most common problem the spacecraft analyst is faced with, however, is the modeling 
of a flexible appendage, such as a boom, antenna, or solar array, which is at a limb end in 
the topological tree model. By making use of equations 107 through 1 10 and equation 
136, the arbitrary body, body X, may first be modeled as a spring-connected rigid body and 
then as a flexible body having several normal modes of vibration. In a step-by-step manner, 
terms which model identical dynamic effects may be equivalenced to each other, and the 
desired equivalence relations, which provide the measure of comparison between tech- 
niques, may be obtained. In the following development, body X is an arbitrary body 
at a limb end. 
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If body X is modeled as a rigid body, and if its 
then, from equations 116 and 129, its relative 
equation 


contiguous body, body J(X), is a rigid body, 
angular velocity is defined by the vector 



m@k-l 


(146) 


The angular velocity of body X, relative to the inertial frame, is given by 



m@k-l 


(147) 


If body X is modeled as a flexible body which is clamped to rigid body J(X), the angular 
velocity of its body-fixed reference frame, relative to the inertial frame, is 

“x = “j ( x) • (148) 

It will be shown that there are three fundamental dynamic characteristics of the body X 
which must be modeled. The modeling technique which best simulates these three charac- 
teristics will provide the most accurate simulation model. The fundamental dynamic 
characteristics are as follows: 

• Relative linear velocity of the center-of-mass of body X 

The linear velocity of the center-of-mass of body X relative to a reference 
frame fixed in the contiguous body, body J(X), must be modeled by both tech- 
niques. The equivalence relation is given by 

^2 k n,xW\,x (149) 

n,\ 

where the left-hand side pertains to a rigid body simulation of body X and the 
right-hand side to a flexible body simulation of the same body X. 

• Relative derivative of the inertia tensor of body X 

The time derivative of the inertia tensor of body X relative to the reference 
frame fixed in body J(X) must be modeled by both techniques. The equivalence 
relation is given by 



n,\ 


(150) 
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• Relative angular momentum of body X about its center-of-mass 

The angular momentum of body X about its own center-of-mass relative to the 
reference frame fixed in body J(X) must be modeled by both techniques. The 
equivalence relation is given by 


$ 





• 

n,\ 


(151) 


To arrive at relations 149 through 151, the system equations of motion, equations 107 
through 1 1 0 and 136, are examined in detail. 

The rotational motion of all nests containing at least one rigid or flexible body is given 
by equation 107. If body X is at a limb end in the nest k-1, the form of the cross- 
coupling term on the left-hand side of the equation is dependent upon how body X 
is modeled. 


On the left-hand side of equation 107, the cross-coupling term which is associated with a 
rigid body model of body X is 


where 



(152) 
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If body X is treated as a flexible body clamped to body J(X), the cross-coupling term on 
the left-hand side of equation 1 07 is 





“n,X 


(t). 


(154) 


The equivalence relation for the cross-coupling term on the left-hand side of equation 136 
is 




*x + G 
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(155) 


where the term proportional to relative angular velocity is carried over to the right-hand 
side of equation 1 36. 
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For small angle relative motion, 
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(156) 
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and the relative angular velocity term 
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Making use of equations 79, 99, 104, and 106, a routine substitution into equation 1 55 yields 
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which readily implies the following relations: 
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The terms on the right-hand side of equation 107, which are dependent upon body deform- 
ation and the inertial angular rates of the bodies contained within the nest k-1 , are given 
by the expression 


Y, [V 5 x + ^x x L x,x + 7 k 


A eS 
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If body X is modeled as a rigid body, then the expression, 

+ 2 -(**»**>)• 


(162) 


k-l 


appearing on the left-hand side of equation 1 07, is carried to the right-hand side of equa- 
tion 136 and must be used for a determination of the equivalence relations. 
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Consider only the contribution of body X contained in the nest k-1. If body X is modeled 
as a rigid body, then it follows from equations 100, 101, 104, 146, and 147, along with 
the fact that body X deformation is zero, that 
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where it is assumed (for simplicity) that there are no momentum wheels in body X. 

If body X is modeled as a flexible body clamped to rigid body J(X), then it follows from 
equations 47, 49, 53, 100, 101, and 148 that 
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The desired equivalence relation is obtained by a substitution of equations 163 through 
169 into equation 161 minus 162 and equivalencing the result to that obtained by a sub- 
stitution of equations 170 through 175 into equation 161 alone. The above substitutions, 
along with a slight rearrangement of terms and deletion of terms quadratic in deformation, 
yield : 
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By equivalencing equivalent dynamic effects in equation 176 on a term-by-term basis and 
applying the vector-tensor identities of the Appendix, the following three equivalence 
equations are obtained: 
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Direct elimination of common factors from both sides of equations 177, 178, and 179 
yields the three fundamental equivalence relations given by equations 149, 150, and 151. 
That is. 
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In an analogous manner, equivalence relations can be derived from the translation equation 
(108) for the nest k-1. The equivalence relation for terms which appear on the left-hand 
side of equation 136 is 
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and the equivalence relation for the terms on the right-hand side of the equation is 
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Direct cancellation of common factors in equations 180 and 181 leads to equations 160 
and 149. That is, 
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where use of the fact that 
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has been made. 
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If body X is modeled as a rigid body, its equation of motion is the rotation equation for 
the nest X-l, equation 107. This equation defines the rate of change of angular momen- 
tum of the body A relative to the inertial origin. 


If body X is modeled as a flexible body clamped to rigid body J(A), its equation of motion is 
the set of simultaneous generalized displacement coordinate equations given by equation 110. 
These equations have the units of translation; however, the generalized displacement coordin- 
ates are used to define both the translational and rotational motion of all elements of the body A. 

Since the units of equation 107 and 110 are different, it is not possible to define equiva- 
lence relations. It should be noted, however, that the coefficient matrix on the left-hand 
side of equation 136 is symmetric for either modeling approach. Hence, the coupling 
coefficients appearing in the body X equation also appear in the equation for each nest 
which contains body X. 


In the right-hand side of the deformation equation (110), a back substitution of original 
definitions provided by equations 48 and 65 yields 
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(184) 


is the force associated with the element i of body X due to its centripetal and Coriolis 
acceleration, relative to an inertial reference, and 
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is the torque associated with the rate of change of relative angular momentum of the ele- 
ment i of body X about its own center-of-mass due to the angular velocity of the body X 
fixed reference frame. 


The other term on the right-hand side of equation 1 1 0, namely 


-m. 




( 186 ) 


is simply the standard spring and viscous damping effects which are normally expected in 
the equations of vibration. 


INTERPRETATION OF RESULTS 

The work presented in this document and in reference 1 provides a complete derivation of 
the equations of motion for a topological tree of coupled rigid bodies, flexible bodies, 
point masses, and symmetric momentum wheels. Furthermore, for the modeling of an 
elastic appendage, equivalence relations have been derived which define exactly how rigid 
body modeling and flexible body modeling procedures purport to describe identical dynam- 
ic effects. 

From an analysis of the derived equations and a knowledge of the digital computation 
which will be required to obtain numerical solutions, several points should be discussed. 
These points are intended to provide the analyst with a set of benchmarks which may be 
used in the development of a simulation model for a particular problem. 


Problem Setup 

The setup of the input data for a coupled rigid body mode] is extremely simple and can 
usually be done in a matter of minutes; furthermore, the input data consists of physically 
realizable quantities which can be readily checked. The setup of the input data for a 
coupled flexible body model is slightly more complex. Several mode-dependent parame- 
ters must be computed which do not have an obvious physical analogy. Depending upon 
the particular problem, these can be time consuming to compute and check. 


Computation Speed 

The internal computation which must be carried out by the computer should increase 
roughly as the square of the number of bodies used. If only first-mode bending vibration 
is to be modeled, the internal computation which will be performed if the body is treated 
as a flexible body is roughly equivalent to that which will be done if it is treated as a 
spring-connected rigid body. 



If higher mode vibrational response is to be modeled, the input of modal data for the flexi- 
ble body model short circuits much of the internal machinery which must be gone through 
every integration step for a coupled rigid body model of the same appendage. In essence, 
a coupled rigid body model of a flexible body is forced to set up and solve an eigenvalue 
problem every integration step, while the flexible body model solves the problem once and 
uses the results to form the required input data. For this situation, flexible body modeling 
is usually mandatory from the standpoint of computation speed. One problem which should 
be recognized is that, for large angle elastic deformation, mode shapes and frequencies based 
upon linear elasticity can be an error source. A coupled rigid body model for this case may 
be necessary in spite of the impact upon computation speed. 

Interpretation of Numerical Solutions 

The interpretation of the results of a coupled body simulation can at times be extremely 
difficult, especially if superfluous degrees of freedom are permitted. The reader should 
be cautioned about using high-degree-of-freedom simulation models. Not only is compu- 
tation speed compromised, but a veil of noise can be created which will mask important 
effects such as slowly growing fundamental instabilities. These remarks apply equally to 
both rigid and flexible body simulation models. 

Choice of Modeling Technique 

The analyst must inevitably decide if the particular problem at hand can be adequately 
simulated as a system of coupled rigid bodies or if some of the bodies must be treated 
as flexible bodies having several natural modes of vibration. Neither modeling approach 
is perfect. 

The rigid body modeling approach permits only fundamental vibrational motion to be 
modeled. This is done by the introduction of a spring constraint which will yield the 
desired natural frequency of vibration for the spring-connected rigid body model. If this 
is acceptable, the equations as derived introduce no further assumptions; they are the 
exact nonlinear equations of motion of the defined coupled rigid body system. 

The flexible body modeling approach permits any number of vibration modes to be used 
to describe the dynamic characteristics of the body. If this approach is used, all terms 
quadratic in appendage deformation are assumed of second order importance and are 
dropped from the system equations of motion. It is emphasized that this assumption 
perpetuates itself throughout the entire set of system equations. It does not, definitely 
not, localize itself in the generalized coordinate equations alone. 

It is up to the analyst to decide whether it is more important for a particular application 
to include higher modes of vibration or to include terms which are nonlinear in appendage 
deformation. The simulation model may be provided with the capability to have one effect 
or the other but not both. 
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The detailed study of the equations of motion which result from both the rigid body and 
the flexible body modeling approach has provided three fundamental equivalence rela- 
tions, namely: 


and 



(149) 


(150) 
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In addition to these relations, the related supplementary relations 


ct 


°x + 


n,\ 


a n,xWA 


n,\ : 


(47) 


and 
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should also be considered. It is essential to realize that the equivalence relations must be 
viewed relative to a common frame of reference. 

To determine which modeling approach is best for a particular application, the analyst 
need only substitute the expected parameters of the problem into equations 149, 150, 
and 151. The side of the equation which yields the most representative description of the 
respective dynamic effects dictates the modeling approach to be used. Furthermore, the 
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right-hand side of the equations provides a guide to the number of mode shapes required 
for representative modeling. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland November 26, 1974 
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APPENDIX 


VECTOR-DYADIC IDENTITIES 

The following vector-tensor identities have been used extensively throughout this docu- 
ment. They are scattered throughout various texts on the subject and collected here for 
the convenience of the reader. 

Let 

"? = vector 

(f = vector 

R = vector 

T = dyad* 

1 = unit dyad 

S’ - skew operator such that for the vector P having components (p t , p 2 , p 3 ), 
relative to a particular frame of reference, 



Making use of the above definitions, the following vector-dyadic identities may be proven: 



PXQ 

= -^(P) • Q 


' -t 

QXR, 

l t 

| - (pXCj) -R 


;sxr) 

\ = q-(rx?) 

P X 

(c?x it) 

ti 

/CH 

i 

PX 

(3xit 

) = (r-p) i-rp • Q 


*Also referred to as dyadic or tensor of rank 2. 
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(?-tJ -r = P • (t-r) 

(txp) • q = t- (?x q) 
p* (qxt) = (?xq) -t 
P • (tXq) = (p-t| XQ 
(pxt) -q = px (t-q) 

(?xt) x q = px (tx q) 
(tX?)xQ = -T (f • Q j + (t • q) P 
PX (^Xt| = - (p • T + Q (p • t) 

Tx(?Xq) = -T- (p^-Qp) 
(px^)x T = -(?Q -Qp) • T 
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